Light Flavour and Heavy Quark Spin Symmetry in Heavy Meson Molecules 
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We propose an effective field theory incorporating light SU(3)-flavour and heavy quark spin 
symmetry to describe charmed meson-antimeson bound states. At lowest order the effective field 
theory entails a remarkable simplification: it only involves contact range interactions among the 
heavy meson and antimeson fields. We show that the isospin violating decays of the A (3872) can 
be used to constraint the interaction between the D and a D* mesons in the isovector channel. As a 
consequence, we can rule out the existence of a isovector partner of the A (3872). If we additionally 
assume that the A(3915) and Y(4140) are D* D* and D* s D* a molecular states, we can determine the 
full spectrum of molecular states with isospin I = 0, | and 1. 

PACS numbers: 03.65.Ge, 13.75.Lb, 14.40.Lb, 14.40.Pq, 14.40.Rt 



I. INTRODUCTION 

The AT (3872) resonance |l[ has opened new perspec- 
tives in hadron spectroscopy. The A" (3872), even though 
it clearly contains a cc pair, does not fit well within 
the standard charmonium spectrum. Disentangling its 
nature requires a more exotic explanation. Among 
the theoretical proposals available, the interpretation of 
the X(3872) as a hadronic molecule [Hill] is especially 
promising and has attracted the attention of the com- 
munity. The motivation behind the molecular picture, in 
which the A" (3872) is a bound state of a charmed meson 
and antimeson, is the striking closeness of this resonance 
to the D°D°* threshold, ra x(3872) = 38 71.68 ± 0.17 Me V, 
to be compared with m D o +m D o, = 3871.84 ±0.20 Me V 
Q. However, the likelihood of the molecular hypothesis 
depends on the J PC quantum numbers of the AT (3872), 
which have not been experimentally determined yet. 
They are either 1 ++ or 2 h [8l-fic| (see also the interest- 
ing analysis of Ref. [HI]), °f which only 1 ++ is compatible 
with a low-lying S-wave bound state. 

The discovery of the AT (3872) has been followed by 
the experimental observation of a series of hidden charm 
resonances above the open-charm threshold, the so- 
called XYZ states. Some of the XYZ states fit well 
within the charmonium spectrum (most notably the 
Z(3940) [l2|), but others do not and may require, just 
like the Y(3872), non-conventional explanations. A few 
might be molecular: several autho rs l l3H l5l| have pro- 
posed the A" (3915) [jj] and Y(4140) fijto be D*D* and 
D* a D* bound states. The F(4260) [3] might even have 
a three body structure (J/* KK) pjjT The Y (4660) 
and AT (4630) J2I have been theorized to be aJVj(980)*' 
molecule [22|,[23| (they may be the same state [23|]). More 
recently, the Belle collaboration has observed two hidden 
bottom resonances, the Z 6 (10610) and Z b (10650) [HlH, 



located just a few MeV away from the BB* and B*B* 
thresholds. Hence, the Z^s may also have a molecular 
nature (26^281 . 

Heavy meson-antimeson molecules are very interest- 
ing objects from the theoretical point of view. As far 
as they are not tightly bound, the meson and antimeson 
will preserve their individuality and will not probe the 
specific details of the short range interaction responsible 
of their binding. Thus, there exists a scale separation 
and molecular states may be amenable to an effective 
field theory (EFT) treatment in which they can be de- 
scribed as mesons interacting via contact interactions and 
pion exchanges (e.g. X-EFT [29], Heavy Meson Molecule 
EFT 30]). In addition, the heavy-light quark content of 
heavy meson molecules implies a high degree of symme- 
try. They are subjected to heavy quark spin symmetry 
(HQSS) [3l1 - l34| . which imposes interesting constraints in 
the heavy meson-antimeson interaction |35j |. As a con- 
sequence, HQSS can be used to predict the existence of 
so far unobserved molecular states [III [28|, [3(| . On the 
other hand, we must not forget the light quark content 
of heavy meson molecules. If we consider q = u,d, s, we 
expect SU(3) flavour symmetry to hold: we can arrange 
molecular states within SU(3) multiplets. In this work 
we will explore the consequences of this symmetry in the 
spectrum of heavy meson molecules. 

The effect of HQSS and SU(3) flavour symmetry 
is to generate relationships among the heavy meson- 
antimesons interactions in different channels. Four pa- 
rameters are enough to describe the twenty-four possible 
S-wave moleculeo That is, we need four data points to 
predict the full molecular spectrum. For this purpose we 
will assume the molecular nature of certain XYZ states 
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1 This is the total dimension, without considering the spin- 
isospin third component multiplicities of the HH space (H = 
D+, D°, D°, D*+,D*° and D*°) or states that are connected by 
a C-parity transformation, such as D a D and D S D that are not 
counted twice. 
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such as the V(3872), X(3915) and F(4140). If the pre- 
dictions turn out to be correct, this will serve as a con- 
firmation that the states used as input are molecular. 
We notice that the A"(3915) could also be a charmonium 
state, as the decay properties are similar to a molecule's, 
with the most important decay channels being J /if? to and 
DD. However, according to Ref. (37|, the charmonium 
hypothesis implies a DD partial decay width of about 
one hundred MeV at least, larger than the experimen- 
tal total decay width (about 30 MeV). In contrast, the 
molecular picture predicts a partial decay width to DD 
of the order of tens of MeV [36[ . We also warn that the 
V(4140) is far from bein g co nfirmed and might not even 
exist: after its discovery nij, it has not been observed in 
subsequent experiments 38, 39]. Finally, we notice that 
SU(2) isospin and SU(3) flavour symmetry are broken 
in the masses of the heavy mesons. From this we can 
derive interesting consequences. In particular, from the 
analysis of the isospin violating branching of the A" (3872) 
into J/if?uj and J /if? p, we will be able to determine the 
strength of the contact interaction in the isovector chan- 
nel and to discard the existence of an isovector or D S D* S 
counterpart of the A (3872) 

The article is structured as follows: in Sect. [TT] we 
briefly review the effective field theory formalism that 
we use for the description of heavy meson molecules. In 
Sect. IHII we consider the isospin violating decay of the 
A (3872) into J /if? p, which can be used to obtain infor- 
mation about the interaction of the DD* in the isovec- 
tor channel. In Sect. IIVI we calculate the location of 
the SU(3)-flavour and HQSS partners of the A (3872), 
A(3915) and V(4140). Finally, we present our conclu- 
sions in Sect. [V] 



II. THE EFT DESCRIPTION AT LOWEST 
ORDER 



A. Overview of the EFT Formalism 

EFTs are generic theoretical descriptions of low en- 
ergy phenomena. They become an adequate tool for 
situations in which a more fundamental description in 
terms of the underlying high-energy dynamics is imprac- 
tical for whatever reasons. Thus they are very useful 
for low-energy hadronic processes where quantum chro- 
modynamics (QCD) is not solvable owing to asymptotic 
freedom and confinement. 

The formulation of an EFT requires the identification 
of the degrees of freedom and symmetries that are im- 
portant for the low energy dynamics. For heavy meson 
molecules, the degrees of freedom are the heavy meson 
and antimeson, plus the pion field. The relevant sym- 
metries, which provide the connection to the underlying 
theory (QCD), are chiral symmetry and HQSS. It is im- 
portant to notice that EFTs are only useful if there ex- 
ists a separation of scales between low and high energy 
physics. If we name the characteristic low energy scale 
at which the EFT is expected to work as Q, and the high 
energy scale as Ao, EFTs allows to construct the ampli- 
tudes as a power series expansion in terms of the small 
parameter xq — Q/^-o 

a=J2 a(u) = Y, xV oA v) , a) 

meaning that if we truncate the expansion at v = v max , 
the a priori error of the calculation will be £g max+1 . The 
ordering principle behind the EFT expansion is power 
counting. It is related to the scaling properties of the 
diagrams contributing to the EFT amplitude A under a 
rescaling of all the Q's 

A^ D HXQ,A )^y D A^(Q,A ), (2) 

where the superscript (D) indicates that we are consid- 
ering a single diagram of order vr> (equivalently, we say 
that D is of order Q VD ). At each order Q v in the EFT 
expansion there is only a finite number of diagrams in- 
volving the low energy fields and their symmetries. 



In this section we review the EFT we use for describ- 
ing heavy meson molecules. The presentation is brief - 
we have already discussed this EFT in previous publica- 
tions [13, HH - and we concentrate on the ideas rather 
than the technicalities. The EFT description involves pi- 
ons and heavy meson/ antimeson fields and the local in- 
teractions among these degrees of freedom that are com- 
patible with the known low energy symmetries, most no- 
tably HQSS and chiral symmetry. A remarkable simpli- 
fication is that pion exchanges are weaker than naively 
expected [3(| and only enter as a perturbation at sublead- 
ing orders. A similar thing happens to coupled channel 
dynamics. Hence, at lowest or leading order (LO), the 
EFT consists on heavy meson and antimeson interacting 
through a contact range potential. 



B. The EFT Potential 

Heavy mesons are non-relativistic and form bound 
states. The first property entails a simplification of the 
EFT description, while the second triggers a series of non- 
trivial changes in the power counting that we will com- 
ment later. The simplification is that the heavy meson- 
antimeson potential is a well-defined quantity. Instead 
of directly expanding the scattering amplitude within 
EFT, we can expand (and truncate) the heavy meson- 
antimeson potential 

v = J2 y (v) + o{x v — +i ) . (3) 

u=0 
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Then we calculate wave functions and observables in the 
standard quantum mechanical fashion. For example, we 
can generate bound states by iterating the EFT potential 
in the Schrodinger / Lippmann-Schwinger equation 0. 

At LO, which corresponds to v = 0, the heavy meson- 
antimeson potential is local 



(p'|F(°)|p)=y(°)(p'-p), 



(4) 



and receives the contribution of two diagrams, a four 
heavy meson vertex and the one pion exchange (OPE) 
potential. The LO potential reads (schematically) 



y(°)(«f) = c( 0) + 



9 2 (a 



q) (b-q) 



<1 ■ 



(5) 



The contact range coupling Cq ^ is a free parameter. Its 
value can be determined from the location of a bound 
state. The OPE contribution depends on the pion decay 
constant ~ 132 MeV and the axial coupling g ~ 0.6, 
which we have particularized for the charmed meson case. 
The sign r\ and the spin operators a and b depend on 
whether the incoming/outgoing heavy meson and an- 
timeson are pseudoscalar or vector (the details can be 
consulted in Ref. (30l|V 

A problematic feature of the EFT potential is that it 
tends to a constant, non-zero value at large exchanged 
momentum. Thus, heavy meson- ant imeson loops con- 
taining the EFT potential will be divergent. We solve 
this issue by renormalizing the EFT calculations. For 
that, we define a regularized potential Va as follows 



(riw> = /(|)G?w) /< 



•A' 



(0) 



where A is an ultraviolet cut-off, p = \p\, p' = \p'\ 
and f{x) a regulator function obeying the conditions (i) 
/(0) = 1, (ii) f(x) — > for x — > oo. For concreteness we 
use the gaussian regulator function 



/Gauss (x) = e 



(7) 



though other choices of the regulator are equally valid. 
The cut-off dependence of the theory is absorbed in 
the counterterms, completing the renormalization pro- 
cess. We comment that there is no universal criterion 
for choosing the cut-off, but as a general rule it does not 
need to be much larger than the high energy scale Ao [H| . 
Here we employ the window A = 0.5 — 1.0 GeV, which is 
of the order of Ap. 



C. Power Counting and Bound States 

The existence of bound states implies changes in the 
power counting of the potential. This can be easily ap- 
preciated by considering the bound state equation 



I*/ 



G (E)V\9 B ), 



(8) 



where \^b) is the wave function, Gq{E) = 1/(E—H ) the 
resolvent operator and V the non-relativistic potential. 
Then it is apparent from the EFT point of view that 
power counting requires the successive iterations of G$V 
to be of the same order 



O(G V) = O(G VG V) . 



(9) 



If we take into account that Gq(E) is of order Q, it 
is clear that the existence of a bound states requires 
that the EFT potential contains a contribution of or- 
der Q -1 . That is, the naive assignment of order Q° to 
the lowest order potential is incompatible with the ex- 
istence of bound states. The solution is to promote at 
least one of the dia gra ms conforming the Q° potential to 
order Q^ 1 [42l. |43|. I45l - l47j . For heavy meson-antimeson 
molecules we usually move the contact range potential to 
Q _1 and left the OPE potential at Q°. The consequences 
of this promotion are that (i) we redefine LO as C _1 and 
(ii) OPE is a subleading order effect, that is, a small 
perturbation over the LO results. As we are only consid- 
ering the LO calculation, we are left with a very simple 
theory that only contains contact interactions. The ex- 
ception is the isoscalar bottom sector, in which OPE is 
stronger than naively expected and is moved to order 
Q^ 1 together with the counterterm (30j . 

To end the discussion, we find it useful to write down 
the eigenvalue equation for the bound states: 



Co (A) 



d 3 q 
(2tt) 3 



2fj, 



A J 



■T 



(10) 



which relates the position of the bound state with the 
strength of the contact interaction Co (i.e. we are as- 
suming that this is the only operator that enters at LO). 
In the equation above f(x) is the regulator function and 
7 2 = —2[iEb, with Eb < the bound state energy and 
fi the reduced mass of the two-body system. Provided 
the cut-off A and the wave number 7, we can determine 
the value of Co- Conversely, from Co and A, we can pre- 
dict the location of a bound state. Further details on the 
eigenvalue equation above and on the wave functions we 
obtain from it can be found in Refs. 13a. l4al. 



2 We notice that this in fact the basis of the Weinberg counting as 
originally formulated in the context of two-nucleon scattering in 
EFT l4Cl |4ll . However, the power counting we obtain is equiv- 
alent to the Kaplan, Savage and Wise proposal \4'A 1-431] . modulo 
the obvious difference that here we are dealing with heavy mesons 
instead of nucleons. 



D. The Heavy Quark Spin Structure of the EFT 
Potential 

According to the previous arguments, at lowest or- 
der in the EFT expansion the potential only contains 
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an energy- and momentum- independent S-wave con- 
tact range interaction (per channel). In a given isospin- 
strangeness sector, naively this translates into six coun- 
terterms, one for each of the possible S-wave configura- 
tions of a heavy meson molecule. However, HQSS reduces 
the number of independent counterterms to two, which 
we call C a and C b ja|. 

Owing to HQSS the LO potential mixes the different 
particle channels. The reason is that the P (P) and P* 
(P*) mesons (antimesons) can be transformed into each 
other by means of a flip of the spin of the heavy quark. 
Nevertheless the total angular momentum and parity of 
a heavy meson molecule is conserved. If we consider the 
following set of basis (with well-defined J p ) 



B(0 + ) = (|PP),|P*P*(0))} , 
B{\+) = {\PP*), \P*P), |P*P*(1))} 
B(2+) = {\P*P*(2))} , 

the LO potential takes the form 

C a vo Cb 
v3 Cb C a — 2Cb y 

C a —C b V%Cb 

"Cb C a V%Cb 

V2Cb V2Cb C a — Cb/ 
C a + Cb ■ 




(11) 

(12) 
(13) 



(14) 



(15) 



(16) 



The notation |p(*)p(*)) is used to indicate a sys- 
tem of a heavy meson/ antimeson containing a heavy 
quark/antiquark, regardless of the light quark content. 
That is, the heavy meson pW is not necessarily the an- 
tiparticle of the heavy antimeson P", as they could con- 
tain different types of light quarks. The number in paren- 
thesis in the |P*P*(J)) states is the total intrinsic spin 
J to which the two vector meson system couples. 

In the specific case in with the heavy meson and an- 
timeson are each other antiparticle, the C-parity is a well- 
defined (and conserved) quantum number. The J p = + 
and 2 + states are now J PC = ++ and 2 ++ states, but 
neither the potentials nor the bases we use to express 
them change 



£(0++) = 6(0+), 
B{2++) =6(2+), 

V W (q,0+ + ) = V hO (q,0 + ), 
V LO {q,2++) = V LO {q,2+). 



(17) 
(18) 

(19) 
(20) 



In contrast, the 1 + states can be further subdivided into 
1 ++ and 1 H states that do not mix under the effect of 
the LO potential. We can write a 1 ++ and 1 H basis 

B(l + -) = j-L (|PP*) + |P*P» , |P*P*(1)}} ,(21) 



(22) 



B(l++) = \ -j= (|PP*) - |P*P» 



for which the LO potential now reads 



C a — Cb 2 Cb 
2Cb C a — Cb 



V LO (q,l + -) 
V LO (q;,l ++ ) = Ca + C 



(23) 
(24) 



As can be seen, the 1 ++ heavy meson- antimeson molecule 
decouples from the two 1 H components. In addition, 
the 1 ++ potential is identical to the 2 ++ one. A 1 ++ 
heavy meson-antimeson molecule implies a 2 ++ HQSS 
partner |36| . 

There is still a significant simplification in the EFT 
potential that we have not discussed yet. In the charm 
sector the EFT description of heavy meson molecules 
is only expected to be valid if the binding energy is 
smaller than a maximum value P ma x of the order of 
150 - 300 MeV (that is, we take P max = with 
7max = Ao ~ 0.5 — 1 GeV and /i = TO£>/2, where tod is 
the D meson mass). This figure is similar to the mass 
gap among the heavy meson-antimeson particle channels: 



A(0+) = M(D*D*) - M(DD) 
A(l+) = M (D*D*) - M(DD*) 



280 MeV, (25) 
140 MeV, (26) 



where A(J P ) denotes the gap in the J p channel. The 
prize of ignoring the couped channel effects is a re- 
duction of the expansion parameter of the theory from 
•Eo = \/B/B m&x to x' = ^/£?/A(J p ), where B is the 
binding energy of the molecular state and we have in- 
cluded the square roots to translate the energy ratios into 
momentum ratios (we are using a non-relativistic EFT). 
As the sizes of P max and A(J P ) are comparable, we can 
safely ignore the coupled channel dynamics without com- 
promising the range of validity of the EFT description of 
the bound states. Moreover, explicit calculations carried 
out in Ref . [HI, |48| for the isoscalar channels confirm that 
the size of the coupled channel effects is the expected one 
in EFT. 

Thus, we can simplify the EFT potential in the + and 
1 + channels to 



V^(q,0+) = C a , 



'pp 
Vp L ,°p,(9,0 



^°p /p p*^ 1+ ) = 



C a ~ 2 Cb , 

+ ^ _ f C a —Cb 
Cb C a 



V£,° p ,(q,l + )=C a -Cb, 



(27) 
(28) 

(29) 

(30) 



where we only need to keep track of the coupled channel 
dynamics in the 1 + p*p/pp* case. If C-parity is a good 
quantum number, the P*P and PP* thresholds coincide 
and the l" 1 component separates from the 1 ++ one: 



LO 



P*P/PP 



p. (9,1+-) =C a -C b 



(31) 
(32) 



V P L 3 /P p,(<7U ++ ) =C a + C b , 
while there is no change in the other channels (except the 
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additional C-parity quantum number) 

V^(q,0 + +) = V^(q,0+), 



'PP 
Fp L , O p,(?>0 H 



'PP 

= v£$.(£o+) 



(33) 
(34) 
(35) 



Now the potentials in the 1+" p*p/pp* and P*P* 
channels coincide. That is, the l" 1 molecular states 
come in pairs, one state per each of the particle chan- 
nels, as happens with the Z;,(10610) and ^(10650) in 
the bottom sector. 



E. The SU(3)-Flavour Structure of the EFT 
Potential 



SU(3)-fiavour symmetry implies that we can organize 
the heavy meson-antimeson states into SU(3)-multiplets. 
As there is a light quark and antiquark (3 and 3 represen- 
tations), heavy meson-antimeson molecules can be orga- 
nized in a singlet and octet representation, 3 <Ei 3 = 1 © 8. 
We can label the heavy meson with a flavour index, that 
is, we write P a , P a , P* and P* a with a = 1, 2, 3, where 
the quark content is 

(Qu\ 

P a ,P: = (Qu, Qd, Qs) and P a , P* a = \Qd\ , (36) 

\Qs) 

and Q = c,b represents the heavy quark. In this notation 
we can construct the singlet and octet representations as 
follows 



\PP,D 



\pp,8;j) 



l 

7! 
i 

71 



PnP a ) 



(^n\PaP b 



(37) 
(38) 



plus the analogous expressions for PP*, P*P and P*P* : 
where the X/s are the Gell-Mann matrices. SU(3) flavour 
symmetry implies that the heavy meson-antimeson inter- 
action distinguishes the singlet and octet representations 

(PP,l\V LO \PP,l) = C a 1] , (39) 
(PP, 8; i\V LO \PP, 8;j) = Cf ) S tJ , (40) 

and do not mix them 



PP,l\V LO \PP,8;j) = 0, 



(41) 



where we have particularized for the case of a PP 
molecule with quantum numbers J PC = ++ and ig- 
nored the particle coupled channel dynamics stemming 
from HQSS. The extension to other particle channels and 
J PC quantum numbers is trivial and only entails the sub- 
stitution of Ca in the equation above by the adequate 
linear combination of Ca and , with fi = 1 or 8. 

Flavour symmetry also conserves the strangeness 5 
and the isospin / of the heavy meson-antimeson system. 
For the singlet representation we have 5 = and 1 = 0. 



For the octet, however, we can have different values of 
the quantum numbers 5 and /, and thus we may find 
more convenient to write them in the basis 



\PP,8;S;I,Mi) 



(42) 



instead of \PP,8;j), where the relation between the two 
basis can be readily obtained (the only technicality in- 
volved is the isospin phase convention for the light anti- 
quarks). In this basis we have 

(PP, 8; S; I; Af 7 | V LO |PP, 8; S'; I'; M'j) 

= C a ^ 5 SS , Sir S Mi m' , (43) 

plus the corresponding expressions for the particle chan- 
nels (PP*,P*P and P*P*). The conservation of S, I 
and Mi is now evident. 

However, flavour symmetry is broken at the level of 
particle masses. In general, this only represents a prob- 
lem for isoscalar states with hidden strangeness (1 = 
and 5 = 0), where the threshold of the particle channel 
that contains the ss light quark pair is considerably above 
the uu and dd channels. For example, in the charm sector 
the energy gap between the D S D S and the D° D° / D + D~ 
thresholds is about 200 MeV, of the order of the maxi- 
mum binding energy we expect to be able to describe 
within the EFT framework. In other cases, the particle 
channels within a given 5 and / channel have approxi- 
mately the same masses. Isospin symmetry breaking is 
of the order of a few MeV and we will not take it into 
account unless the molecular state we are describing is 
very shallow, as in the X(3872), a case we will consider 
in detail at the end of this section. 

For the treatment of SU(3) breaking we define 



|PP [°] 



1 

71 



(Ipp 1 ) + |p 2 p 2 )) , 



|p s p s ) = |p 3 p 3 ) 



(44) 
(45) 

that is, we distinguish between the SU(2) isoscalar (I = 
0) state and the hidden strangeness state, and the indices 
1, 2 and 3 refer to the representation of Eq. (|3"6")l . In 
this basis we can rewrite the singlet and octet isoscalar 
representations as 

|PP,1) = -4 \V2\PP[0]) + \P S P S ) 



V3 

|PP,8;0;0,0) = i 



|PP[0]) - V2\P S P S 



(46) 
(47) 



By inverting these relations, we find that in the basis 

6={|PP[0]),|P S P S )} , (48) 
the flavour symmetric interaction can be written as 



yLO = 



2 r W 
3 U 



3 



1^(8) 

- c^' 



V2 
3 



3° Q 



(i) 



2 r (8)" 



,(49) 



where flavour symmetry is broken owing to the energy 
gap among the |PP[0]) and |P S P S ) thresholds. As com- 
mented in the previous paragraph, the energy gap is large 
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enough as to justify the treatment of the isoscalar and 
hidden strangeness channels as uncoupled. 

In terms of notation we find that working in the SU(2) 
isospin basis is easier than in the SU(3) one. If we define 
the SU(2) isovector state (with Mj = 0) as 



\PP[1] 



1 

71 



(IAP 1 ) - |P 2 F 2 )) , 



(50) 



we can rewrite the EFT potential in the SU(2) isospin 
basis as 



(PP[0}\V LO \PP[0}) = C 0a . 
(PP[l]\V LO \PP[l]) = C la . 



(51) 
(52) 



plus the corresponding expressions for the other particle 
(PP*,P*P and p*p*) j spin / C-parity combinations, 
which require the inclusion of the new counterterms Cob 
and C\b- 

The counterterms Co a arid C\ a introduced in Eqs. (|51|) 
and (|52[) are related to the ones in the original SU(3)- 
multiplet basis by 



c 0a = \cp + \cw, 
c la = cp. 



(53) 
(54) 



In the SU(2) basis, the interaction in the \P S P S ) channel 
reads 



(P.P.W 



LO| 



1 



(Coa + Ci a ) 



(55) 



that is, the average of the isoscalar and isovector contact 
range potentials. 

Finally, we consider the X(3872), where isospin break- 
ing is important as this bound state is especially shal- 
low. The energy gap between the D°D*° and D + D*~ is 
8MeV, which is smaller than the binding energy of the 
X(3872) (about 4MeV in the isospin symmetric limit, 
or, equivalently, almost at the D°D*° threshold). What 
we do then is to treat the neutral (D°D*°) and charged 
(D + D*~) channels independently, which are related to 
the isoscalar and isovector channels by 



\DD* [0]) = i [\D°D*°) 
v2 

\DD* [1]) = -L [\D°D°*] 
v2 



\D + D* 



\D+D* 



(56) 
(57) 



and assume that isospin symmetry is only broken at the 
level of the masses. If we consider the basis of physical 
states 



B 



X(3872) 



1 

71 



D°D*°) - \D*°D )) , 



V2 



(\D+D*-) - \D*+D-)) 



we find that we can express the EFT potential as 



V LO = - 



Co + Ci Co — Ci 
Co — C\ Co + Ci 



(58) 



(59) 



where Co and C\ stand for 

Co = C 0a + C b , C\ — Ci a + C\b , (60) 

that is, the linear combination of isoscalar and isovector 
counterterms corresponding to the 1 ++ channel. 



III. ISOSPIN SYMMETRY VIOLATION IN THE 

X(3872) 

The Belle collaboration reported for the first time the 
decays of the X(3872) into the (isoscalar) J/^>tt + tt^tt° 
and the (isovector) J/tyn+n^ channels The non- 
negligible size of the latter one hints to the existence 
of isospin breaking terms in the dynamics that govern 
these decays. The size of the violation is remarkable, as 
indicated by the branching ratio [49j 



B x = 



T(X(3872) -> J/^TT+TT-TT ) 
r(X(3872) -> J/tf 7T+7T-) 



= 0.8 ±0.3. 



(61) 



where the central value is even smaller than one. The 
most natural explanation for the large ratio is the isospin 
breaking generated by the mass difference of the neu- 
tral (D°D *) and charged (D + D*~) channels in the 
X(3872) [5(3, [5l|, which would not have a definite isospin. 
In this picture, at short DD* distances, the X(3872) 
would be a linear combination of I = and 1 = 1 com- 
ponents. The J/^3tt and J/\l/27r decays would proceed 
through the isospin invariant coupling of these final states 
to the appropriated short distance X(3872) isospin wave 
function components [5lj . 

We expect the 2tt and 37r decays to happen via an 
intermediate p and u) meson. In this regard it is interest- 
ing to notice the careful analysis of Hanhart et al. [ll| , in 
which the branching ratio Bx is translated into the more 
convenient ratio 

Rx M(X^Jj¥u) °- 26 - 005 ' (62) 



that, instead of the decay widths, involves the ampli- 
tude of the X(3872) to decay into J/^ and the p or ui 
mesons. In Ref. [ll|, the transition amplitudes for the 
decays of the X(3872) into J/\& and a vector meson are 
parametrized as 



M(X^ J/^V)=g x{ v)fx{p), 



(63) 



where V — w, p and fx is a Blatt-Weisskopf barrier factor 
that depends on the momentum p of the J/\P in the X 
rest frame. For the 1 ++ assignment, we have fx(p) = 1 
yielding the value quoted in Eq. (R)2"j) for Rx ■ Other J PC 
assignments of the X(3872), in particular J PC = 2~+, 
will imply a different form for the decay amplitude and 
consequently a change in the value of Rx (see Ref. [ll| 
for further details). 

If we assume that the X(3872) is a DD* molecule with 
quantum numbers 1 ++ , we can calculate gx(u) and gx(p) 
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in terms of the analogous g u and g p couplings for the free 
DD* meson pair, which are given by 

M (DD* (IS) -> J/* u) = g u , (64) 
M{DD*(IV)^ J/y P ) = g p , (65) 

where IS and IV indicate the isoscalar and isovector 
configurations respectively. When the mesons are bound 
in the A (3872), we can dehne the isoscalar and isovector 
wave functions as the brackets 

(X\DD*(IS)q)=y x{IS) (q), (66) 
(X\DD*(IV)q) = y X (iv)(q) , (67) 

where we have explicitly indicated the relative momen- 
tum q of the meson-antimeson pair. Now, we can re- 
late the g^ and g p couplings with the gx(u) an d 9x(p) 
ones by noticing that we can rewrite the decay ampli- 
tude M(A —> B) as the bracket (B\A). Then, we insert 
the two-body identities 



lis 



rf 3 



(2tt) 3 
d 3 q 
(2^T 



DD* (IS) q) {DD* (IS) q\ , (68) 
3l DD*(IV)q)(DD*(IV)q\, (69) 



for the isoscalar and isovector channels within the brack- 
ets and end up with 

M (X -> J/* u) = g X (u) = 9u *x(/S) , (70) 
M(X^ J/^p) = g x{p) =g P ^ X (iv), (71) 

where 'i'xtis) an d ^x(iv) are defined as 



X(IS) 



X(IV) 



^1 

(27T) 

(2tt)* 



3 *X(/S) 



*X(7V)(?) , 



(72) 
(73) 



that is, they are the coordinate space wave functions at 
the origin (more properly, around the origin, as the regu- 
lator smears the wave functions) 0. Putting all the pieces 
together, we find that the ratio Rx can be expressed as 



Rx 



9x{ P ) _ gp_ ^x(iv) 
9x(w) 9u 4f X (is) 



(74) 



However, we are not working in the isospin basis, but 
rather in the particle basis. If we denote the neutral and 
charged channels by the subscripts and 1, we have that 

*xo = (^X(IS) + ^x(iv)) , (75) 

*xi = 77= (^x(is) - ^x(iv)) , (76) 



from which we can rewrite Rx as 



Rx = 



g P ^xo - vfrjri 



9uj * 



xo 



(77) 



xi 



At this point we notice that if we know Rx , the binding 
energy of the A (3872) and the g p j g^ ratio, we can deter- 
mine the contact range potential that binds the A (3872) 
(which in turn determines whether the A (3872) has an 
isovector partner) . Of these quantities the only unknown 
is the g p j 'g u ratio, which can be determined from the 
SU(3) relation 



-y/2i 



9 P ~ 9lo 

and the OZI rule, that implies 



9p,9u > 94> > 



(78) 



(79) 



which amounts to ignoring the decay of the A (3872) into 
J/^(f>. Thus, we can conclude that gp/g^ — 1 and that 
the Rx ratio simplifies to 



Rx 



$>XQ 



XI 



xo 



(80) 



A'l 



The next step is to solve the bound state equation for 
the A (3872). We consider that the A (3872) is a DD* 
molecule with quantum numbers J PC = 1 ++ , where we 
can distinguish between the neutral (D°D*°) and charged 
(D + D*~) components of the wave function. That is, 
there are two channels in the bound state equation. We 
regularize the LO potential of Eq. ([BUI with a Gaussian 
regulator function and take a cut-off A = 0.5 — 1.0 GeV. 
The couplings Co and C\ (or equivalently Co a + Cob and 
Cia+Cib) are determined by the condition of reproducing 
the location of the A(3872) state and the p/oj ratio Rx, 
from which we obtain the values 



C = -1.693±S;8Uj fm 2 



-O.OSt^fm 2 



(-0.731™ fin-), (81) 
(-0.373™ fm 2 ), (82) 



for A = 0.5 GeV (lGeV). The couplings indicate that 
the strength of the interaction in the isovector channel 
is weaker than in the isoscalar one. In particular, the 
isovector coupling C\ is not strong enough to generate a 
second bound state, the (mostly) isovector partner of the 
A(3872). 

We hnd it worth commenting the comparison of our re- 
sults for the counterterms with the hidden gauge model 
of Gamermann et al. [13, [EH , in which the value of the 
counterterms are determined from meson exchange sat- 
uration. While the hidden gauge model predicts C\ = 
(as the contribution from the p and uj mesons cancel 
out) and Rx — 0.136 (for a sharp cut-off regulator and 
A = 653 MeV, see Ref. [H| for detail^ ), we obtain a 



4 Note that, when C\ = 0, once the regulator/renormalization 
3 Alternatively, we can employ the coupling versus wave-function procedure and the cut-off is fixed, Rx is totally determined by 



language of Ref. 15 ill to arrive at the same result. the mass of the X(3872). 



8 



small (yet important) contribution to C\. As can be 
seen, this small contribution is necessary to fine tune the 
isospin violating branching ratio of the V(3872) decays 
to its exact value. 



IV. THE SU(3) AND HQSS PARTNERS OF THE 

X(3872) 

If we are able to determine the value of the countert- 
erms of the LO EFT, we can calculate the location of 
the molecular partners of the X(3872). There are four 
parameters. Isoscalar molecules (I = 0, S = 0) without 
hidden strangeness can be described with Co a and Cob. 
Isospinor (I = 1/2, S = ±1) and isovector (7 = 1, S = 0) 
states are in turn determined by C\ a and C\ . Finally, 
for molecular states with hidden strangeness, the con- 
tact range interactions are the average of the isoscalar 
and isovector ones. 

We fix two of the counterterms from the location of 
the X(3872) resonance and its isospin breaking branch- 
ing ratio, as explained in the previous section. The re- 
maining two require the identification of two partners of 
the V(3872)._We have chosen the V(3915) El as a 0++ 
isoscalar D*D* molecule and the Y(4140) [13] as a ++ 
D* S D* S molecule, guided by its apparently dominant de- 
cay into J/^cj). We notice that these identifications were 
proposed for the first time in Refs. [13l - ll5| . 

The ++ assignment for the quantum numbers of the 
V(3915) and Y(4140) deserve some discussion. On the 
one hand, the 2 ++ option is excluded: HQSS fixes the lo- 
cation of the isoscalar 2++ D*D* partner of the V(3872) 
in the vicinity of 4012 MeV, far away from the 3915 MeV 
region. If we additionally consider the isospin breaking 
decays of the V(3872), we can determine that the isovec- 
tor 2++ D*D* and the isoscalar 2++ D* S D* S molecules 
do not exist. On the other the 1 H option can also be 
discarded from the decays of the V(3915) and F(4140) 
into J/ty uj and J/ty 4>, requiring a positive C-parity state. 
Thus we are only left with J PC = ++ . 

However, the choice of the X(3915) and Y(4140) states 
as input is not entirely free of problems. The first is the 
binding energy of the A(3915), approximately 100 MeV. 
Around this binding there is power counting transition 
- one pion exchange changes from perturbative to non- 
perturbative - though the critical binding energy that 
marks this transition is not known exactly and can hap- 
pen at larger bindings than expected [3(| ■ In this regard 
the investigations of Ref. [36[ shows (by performing the 
explicit calculations) that we are still in the perturba- 
tive regime for the A(3915) and we can employ a con- 
tact theory to describe it at LO. The V(4140), with a 
binding energy of 80 MeV, is not affected by this issue 
because the D*D* system cannot exchange a single pion. 
The second problem is the debatable experimental sta- 
tus of the F(4140). After its first and only observation 
by the CDF Collaboration |17|. subsequent experiments 
have failed to find it [H, l39j] . This means that the con- 



sequences derived from the assumption that the V(4140) 
exists and that it has a molecular structure should be 
taken with a grain of salt. 

Now we fix the counterterms to the location of the 
V(3915) and V(4140). For convenience we define first 
the linear combinations of counterterms: 



C 2 

c 3 



Coa — 2Cob 

1 



(83) 



((Coa - 2C ob ) + (Ci„ - 2C lb )) , (84) 



where C 2 / C3 is a convenient way to write the LO po- 
tential of the X(3915) / Y(4140) channel. We obtain 



C 2 = -6.710 fm 2 



C, 



-5.915 far 



(-1.611 fm 2 ), (85) 
(-1.459 fm 2 ), (86) 



for a gaussian regulator A = 0.5 GeV (lGeV). We can 
transform the values of the combinations above (and the 
corresponding ones for Co and Ci ) to the standard coun- 
terterm representation, yielding 



Coa = — 
Cob = + 
Cla = 

(hh = +L<ks 



3.3661°;°* fm 2 



(87) 



1 673 +0 012 fm 2 

1 -°'°-0.008 1111 



'•76™fm 2 



1.02 4 ™fni 2 ) 
(+0.2931^2 fm 2 ) 
(-0.684^ fm 2 ), (89) 
(+0.31i;°°33fm 2 ),(90) 



for A = 0.5 GeV (lGeV), where the error comes from 
the uncertainty in Rx (a negligible effect in the isoscalar 
channels, but important in the isovector ones). 

We notice that there is an additional error source: the 
violations of HQSS due to the finite charm quark mass. 
This effect can be taken into account by noticing that the 
EFT potential has a relative uncertainty of the order of 



V, 



LO 

(tnq=tn c ) 



= V, 



LO 



(1±^2£), 



(91) 



with respect to the exact heavy quark limit. Taking a 
value of around 1.5 GeV for the charm quark mass and 
Aqcd ~ 200 MeV, we should expect a 15% violation 
of HQSS for the LO contact range potentials. We no- 
tice that the two error sources are independent: we can 
compute the total error by adding the partial errors in 
quadratures. 



A. The SU(2) Isoscalar (J = 0) Partners 

We begin with the SU(2) isoscalar sector, in which we 
ignore the hidden strange components. The states are 
determined by the counterterms Co a and Co&. We do not 
take into account particle coupled channel effects as they 
are subleading, as explicitly checked in Ref. [36j . There is 
one exception though, the 1 ++ and 2 ++ channels, where 
the mass gap between the neutral and charged channels 
(8 and 6 MeV in each case) is similar in size to the bind- 
ing energy in the isospin symmetric limit (4 and 5 MeV) . 
This suggests that we may treat the neutral and charged 
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channels as explicit degrees of freedom of the theory. We 
note that the inclusion of isospin violation is the only dif- 
ference with respect to the previous analysis of Ref . (3(| . 
The spectrum of molecular states is presented in Table fl] 
As can be seen, isospin violation is a small perturbation 
over the former predictions of Ref. (3(| (though it is still 
crucial to describe isospin violating decays properly). 

We, however, would like to make a few remarks. The 
first one concerns the 2 ++ state. The central values of the 
counterterms as given in Eqs. (j8"7) - (|5fJl) predict that the 
2 ++ state lies very close (less than 1 MeV) to the lowest 
energy threshold, i.e. D*°D*°. When we decrease the 
strength of the potential to account for the uncertainties 
of our approach, the pole reaches the neutral threshold 
and then bounces back into the second Riemann sheet. 
That is, the state becomes virtual (instead of bound). In 
any case, the existence of the pole will strongly influence 
the amplitude at threshold. 

The second comment is about the ++ DD channel. 
As can be seen in the Table HI this state is bound by 
about 20-25 MeV. For simplicity, we have used in this 
channel the isospin symmetric limit. Yet the D°D° — 
D + D~ threshold gap is around 9 MeV, and it might make 
sense the explicit consideration of isospin breaking in this 
channel. However, as in the V(3872) and V(4012) cases, 
the effect is rather small, justifying the validity of the 
isospin symmetric limit for the spectroscopy problem. 

Finally, there is a remark concerning the possible ef- 
fect of hidden strange channels in the dynamics of the 
V(3915) state. As we will see, there is a ++ D S D S state 
at 3925 MeV (see Table HV|) that in principle can mix with 
the V(3915). However, explicit calculations show that 
the influence of the hidden strange channel is numeri- 
cally marginal: it moves the position of the predicted 
molecular states by a small fraction of a MeV (usually 
AM ~ 0.1 - 0.2 MeV), a tiny effect compared to other 
error sources. The reason lies in the transition potential 
between the isoscalar ++ D*D* and D S D S channels: 



(D*D*\V hO (0+ + )\D s D s ) = ] j~(C 0b -C lb ), (92) 

which turns out to be numerically quite small, since we 
accidentally find C ob ^ C u , see Eqs. (|BTl) - (j9"0")) . 



B. The Isospinor (I = \) Partners 

The isospinor molecules are different in the sense that 
they have not well-defined C-parity, as they are not 
bound states of a heavy meson and its antimeson. In 
general this poses no problem (the formalism is identical 
to the one in the previous case) except for the 1 + D S D* 
and DD* molecules. The D S D* and DD* thresholds are 
separated by only 2 MeV and require a coupled channel 
treatment. The LO potential reads 



where channel 1 (2) corresponds to D S D* (D*D). We 
find that this coupled channel potential only leads to one 
bound state (of a maximum of two). We list this and 
the other isospinor molecules in Table UH where we have 
considered only the strangeness one states. The spectrum 
is identical for the strangeness minus one sector. 

There is a total of four states, one of them - the + 
D S D molecule - almost at threshold. Actually, there ex- 
ists a small violation of the third component of isospin, 
because of the different masses of the D~ (D*~) and D° 
(D*°) mesons. In general the isospin splitting is smaller 
than the other errors in the calculations and we ignore it. 
The exception is the + D S D molecule: in Table |TT] we 
only report the I3 = —1/2 component of the molecule, 
which is being formed by the DfD~ interaction 0. Simi- 
larly to what happens with the isoscalar 2 ++ state, with 
the central values of the counterterms in Eq. (|87|) - (r90l) we 
predict a + state that lies very close to threshold. When 
we decrease the strength of the potential, the molecule 
approaches the threshold and finally becomes a virtual 
state, at least within our simple scheme 0. In the real 
world, it might well happen that in this situation the 
state could become a narrow resonance, placed very close 
to threshold, and decaying into the D+D~ pair instead 
of a becoming a virtual state. In any case, we predict 
the existence of some structure close to threshold. The 
upper errors for the mass of this state, as quoted in Table 
HI! just account for the distance between its central mass 
value and the DfD~ threshold. 



C. The Isovector (/ = 1) Partners 

The spectrum of the isovector molecules is similar to 
the isospinor one. The reason is that the LO potential 
is identical in both cases, with the exception of the 1 ++ 
and 2 ++ molecules owing to isospin violation, as we com- 
mented previously. We find four molecular states of a 
possible total of six that we list in Table Hill The location 
of the states is similar to the isoscalar sector, only that 
they are a bit less bound. The two missing states would 
correspond to the isovector partners of the V(3872) and 
X(4012). 

We also mention that there is a complication with the 
isovector ++ DD molecule, which lies very close to the 
DD threshold. This implies that explicit isospin breaking 
should be taken into account. However, we notice that 
this is only necessary for the I3 = component of the 



V 



LO 



Cla —Clb 
-Clb Ci a , 



(93) 



The dynamics of the other state member of the isospin doublet 
is similar, being its mass just shifted by about 5 MeV. This shift 
is due to the mass difference between the D~ and the D~£ D° 
pairs. 

A resonance state in one channel is usually associated to a barrier 
in coordinate space which is not reproduced by a constant (in 
energy) potential. The situation may be different when we have 
coupled channels , in which case one of the channels can decay 
into the other. 
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Hfi 


2S+l Lj 


V c 


£ (A = 0.5 GcV) 


E (A = 1 GeV) 


Exp [7] 


Threshold [MeV] 


0++ 


DD 


'So 


Coa 


3709 ± 10 




- 


3734.5* 


1++ 


D*D 




Eqs. (55JI and (|BDJ 


Input 


Input 


3871.6 


3871.8/3879.9 


1+- 


D*D 


3 s 1 


Coa — Cob 


3815 ± 17 


38211^3 




3875.9* 


0++ 


D*D* 


'So 


Coa ~ 2 Co6 


Input 


Input 


3917 


4017.3* 


1+- 


D*D* 


S S 1 


Coa — Cob 


3955 ± 17 


3958^27 


3942 


4017.3* 


2++ 


D*D* 


5 S 2 


Eqs. dSSJl and {60]) 


4013!"^ 


40131^ 




4014.0/4020.6 



TABLE I. Predicted masses (in MeV) of the SU(2) isoscalar HQSS partners of the X(3872) resonance for two different values 
of the Gaussian cutoff. We use as input 3871.6 MeV, 3917.4 MeV and 4140 MeV for the V(3872), X(3915) and V(4140), 
respectively. In the last column, we also give the thresholds for the different charge combination channels (P° P° / P + P~) or 
the threshold value that results when we neglect isospin breaking effects and we employ charge averaged masses. In this latter 
case, we mark it with an asterisk (*). Errors in the predicted masses are obtained by adding in quadratures the uncertainties 
stemming from the two sources of systematic errors discussed at the end in Subsect. IIVI errors in our determination of the LO 
counterterms in Eqs. (|87[) - (|90|) an d variations obtained by modifying the strength of the contact interaction in each channel 
by ±15%, which corresponds to the expected violations of HQSS for the charm quark mass, ft : the upper error of the 2 ++ 
state mass deserves a detailed discussion that can be found in the text (Subsect. IIV A|) . 





HH 


2S+l Lj 


Vc 


E(A= 0.5 GeV) 


E (A = 1 GeV) 


Exp [7] 


Threshold [MeV] 


0+ 


D+D- 


'So 


Cla 


3835.8t3'j 


3837.71^3 




3838.1 


1+ 


D B D*,D*D 


3 S 1 


Eq. IJ93J 


3949 ± 13 


3957±^ 




3977.15 t ,3979.55 t 


0+ 


D* S D* 


'So 


Cla ~ 2 Clb 


4056 ± 22 


4061±» 




4120. 9 1 " 


1+ 


D* S D* 


3 s, 


Cla — Clb 


4091±il 


4097±2o 




4120. 9 1 " 


2+ 


D*D* 


5 s 2 


Cla + Clb 











TABLE II. Predicted masses (in MeV) of the isospinor (J = |) HQSS partners of the V(3872) resonance, for two different 
values of the Gaussian cutoff. The meaning of the quoted errors in the table is the same as in Table U We also give the relevant 
thresholds (in MeV) for each channel. We use isospin third component averaged masses for those cases marked with a f symbol. 



isospin triplet, where the neutral (D°D°) and charged 
(D + D~) channels are to be found. If we consider in- 
stead the ^3 = ±1 states (corresponding to the particle 
channels D + D° and D°D~), there is no mixing with the 
isoscalar components and no need for coupled channel 
dynamics. Thus, for simplicity, we have decided to re- 
port only the I3 = +1 state in Table Unl 



D. The Hidden Strange Partners 

Finally, we consider the molecules with hidden 
strangeness. In this sector the strength of the LO po- 
tential is the arithmetic mean of the isoscalar and isovec- 
tor one. This means that if there is a bound state in 
the isoscalar and isovector sector, we can be confident 
about the existence of a heavy meson molecule contain- 
ing a ss quark-antiquark pair. Conversely, if we consider 
the A"(3872) and A"(4012) molecules, the fact that they 
have no isovector partners is a strong hint that there will 
be no hidden strange partner either. The four states we 
obtain are listed in Table HVl 



V. CONCLUSIONS 

In this work we have shown how the heavy quark 
spin and light SU(3) flavour symmetries constrain the 
charmed meson-antimeson interaction (HH — ¥ HH, be- 
ing H = D + , D° 7 D®, D*+,D*° and Df). This has been 
done within the EFT framework, where the heavy meson 
interactions can be easily arranged from more to less rel- 
evant thanks to power counting, the ordering principle 
behind EFT. The bottom line of the EFT approach is 
that contact interactions (i.e. four meson vertices) domi- 
nate the low energy dynamics of heavy meson molecules. 
In turn, SU(3) flavour symmetry and HQSS reduce the 
number of contact interactions from twenty four to only 
four. The approach we advocate is actually an extension 
of the ideas of Ref. 36] for the SU(2) isoscalar channel 
to the isospinor, isovector and hidden strange sectors. 

The four counterterms can be determined as follows: 
we fix three of them by identifying the AT(3915), V(4140) 
and X(3872) resonances as molecular states. In partic- 
ular we consider the X(3915) and V(4140) to be a ++ 
isoscalar D*D* and D* S D* S molecules respectively, while 
the AT (3872) is a 1 ++ isospin admixture DD* state. Fi- 
nally, the fourth counterterm is determined from the 
isospin breaking branching ratio of the X(3872) into 
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HH 


2S+l Lj 


Vc 


£ (A = 0.5 GeV) 


E (A = 1 GeV) 


Exp [7] 


Threshold [MeV] 


0++ 


D+D° 


'So 


Cla 


3732.5t^ 


3734.31^6 


- 


3734.5 


1++ 


D*D 


3 S1 


Eqs. ||SSJ and JHSJ) 


- 


- 


- 




1+- 


D*D 


3 Sl 


Cla — Clfc 


3848tii 


3857^8 




3875.9* 


0++ 


D*_D* 




Cla — 2 Clb 


3953 ± 22 


3960lg 




4017.3* 


1+- 


D*D* 


3 S1 


Cla — Clb 


3988 ± 13 


39951^ 




4017.3* 


2++ 


D*D* 


5 & 


Eqs. ((59j) and (J60j) 











TABLE III. Predicted masses (in MeV) of the SU(2) isovector HQSS partners of the X(3872) resonance for two different values 
of the Gaussian cutoff. The meaning of the quoted errors in the table is the same as in Table [U We also give the relevant 
thresholds (in MeV) for each channel. For the D* D and D*D* cases, we give the threshold that would correspond to the zero 
isospin third component and calculated with charge averaged masses. 





HH 


2S+l Lj 


V c 


E (A = 0.5 GeV) 


E (A = 1 GeV) 


Exp [7] 


Threshold [MeV] 


0++ 


D S D S 


'So 


\ (Coa + Cla) 


3924±? 


3928±ii 




3937.0 


1++ 


D* S D S 


3 Si 


\{Coa + Cla + Cob + Clb) 










1+- 


D* S D S 


3 Si 


Tj(C()a + Cla — Cf)b — Clb) 


4035 ± 15 


4040t24 




4080.8 


0++ 


DID*. 


'So 


\ (Coa + Cla — 2 C()b — 2 Clb) 


Input 


Input 


4140 


4224.6 


1+- 


D*D* S 


3 Si 


2"(Coa + Cla — Cob — Clb) 


4177 ± 16 


4180±24 




4224.6 


2++ 


D*D* S 


5 s 2 


\(p0a + Cla + C()b + Clb) 











TABLE IV. Predicted masses (in MeV) of the hidden strange isoscalar HQSS partners of the X(3872) resonance for two 
different values of the Gaussian cutoff. The meaning of the quoted errors is the same as in Table [T] We also give the relevant 
thresholds (in MeV) for each channel. 

interaction vanishes. Here, we improve on that and take 
advantage of the experimental ratio to constraint the 
non- vanishing interaction in the isovector channel. 

Once we have fixed the counterterms, we have estab- 
lished the existence and the location of up to a total of 15 
molecular partners of the A(3915), V(4140) and A(3872) 
states, see Tables HllIVl These predictions are subjected 
to a series of uncertainties, in particular the approximate 
nature of HQSS (especially in the charm sector). We have 
estimated the size of these corrections and concluded that 
the HQSS pattern of molecular states is stable, though 
the exact location of the states can change by a few tens 
of MeV in certain cases. 

Actually, the family of D^D^*> states we theorize de- 
pends on the assumption that the A'(3872), A(3915) 
and y(4140) resonances are molecular. While in the 
A (3872) case the molecular interpretation is compelling 
and widely accepted, the A(3915) and V(4140) states are 
merely compatible with it. Regarding the A(3915), it is 
interesting to notice that the size of the decay width of 
this resonance is difficult to conciliate with the hypoth- 
esis that it is a charmonium [371 ]. This observation en- 
hances the prospects that the A (3915) may be a molecule 
after all. However, with the V(4140) we have a more 
serious problem: this resonance is far from being con- 
firmed experimentally. Thus not all the states we predict 
are equally likely. Predictions derived from the A (3872) 
should be regarded as more solid than those depending 
on the A (39 15), which in turn are less speculative that 



J/*2tt and J/^3tt. We notice that the A(3872) is 
weakly bound (~ 0.2 MeV) and lies extremely close to 
the D°D°* threshold. Its binding energy is much smaller 
than the mass splitting between the charged and neu- 
tral components of the A (3872) (<~ 8 MeV). Hence, we 
have taken into account these degrees of freedom explic- 
itly. We work within a scheme where all interactions are 
isospin invariant, but where the A(3872) does not have 
a well defined isospin as a consequence of the mass 
and kinetic terms of the DD* Hamiltonian [5(J [H| • The 
isospin breaking of the masses naturally explains the ap- 
pearance of an isospin violation (and its quantity) in the 
decays of the A (3872), as pointed out in Refs. [50|, H3 
for the first time. However, in the aforementioned works 
the experimental branching ratio (49j 

r(A(3872) -> J/tt.W) 
r(A(3872) -> J/^tt+tt-) 

is not perfectly reproduced @. The reason is that the 
authors of Refs. [50, uLUl assume that the isovector DD* 



7 The isospin mixing depends on the relative distance between the 
pseudoscalar and vector mesons. 

8 Notice, however, that Refs. [50t l5ll used an older experimental 
determination of the branching ratio (1.0 ± 0.3 ± 0.4) with much 
larger errors. It is only with the updated branching ratio that 
their explanation fails at the fine quantitative level. 
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the ones obtained from the Y(4140). In this regard, as 
stressed in Ref. |36j, the 2 ++ D*D* isoscalar partner of 
the X(3872) is still the most reliable prediction of the 
present work, followed by the other isoscalar states. If 
in the future we count with clearer molecular candidates 
than the A(3915) and, particularly, the Y(4140) reso- 
nances, they could be included in the current scheme 
instead of the later ones, helping to achieve more robust 
predictions. Conversely, the observation of any of the 
states predicted here can serve as proof of the molecular 
nature of the previous mentioned resonances. 

Finally, we find it interesting to compare our results 
with those of the hidden gauge formalism, another the- 
oretical approach for the study of hidden charm states. 
While the spectrum of the isoscalar molecules in the hid- 
den gauge is similar to the one we obtain (with the no- 
table exception of the 2 ++ state), fewer poles are re- 
ported in the I = 1/2, I = 1 and hidden strange 
sectors [53T4551 ] . If we consider the case of two heavy 
pseudoscalar mesons, Gamermann et al. [HH predict an 
isoscalar ++ DD state in the vicinity of 3700 MeV. We 
can identify this state with the AT(3710) DD molecule 
we obtain in the present work. However, they do not 
predict the existence of the isospinor, isovector and hid- 
den strangeness partners of the A (3710). The reason is 
that in Ref. [53[ the interaction in the isovector chan- 
nel is zero. The same comments apply to the extension 
of the hidden gauge formalism to pseudoscalar-vector 
molecules [54j, where it is found a counterpart of the l" 1 
DD* isoscalar molecule that we obtain at 3820 MeV at a 
slightly higher energy (3840 MeV) but no I = 1/2, I = 1 



or hidden strangeness states. Last, in the case of two 
heavy vector mesons the hidden gauge predicts a series of 
isoscalar ++ , l" 1 and 2 ++ D*D* states, plus a few non- 
isoscalar ones (55[. In the isoscalar sector, the ++ and 
l" 1 D*D + resonances are located in the region around 
3943 and 3945 MeV respectively, not very different to the 
masses we use (3917 and 3955 MeV). However, there is a 
striking difference in the mass of the 2 ++ D*D* isoscalar 
state, which in Ref. [55] happens at 3922 MeV (instead of 
4013 MeV). The reason for this value, which is incompat- 
ible with HQSS, is the remarkably strong vector-vector 
interaction that is obtained in the hidden gauge model. 
Probably as a result of this strong interaction, Ref. [55[ 
also report the existence of broad 2 ++ isovector D*D* 
and hidden strange D*D* resonances, with widths above 
100 MeV, and masses of around 3910 and 4160 MeV, 
respectively. These states are difficult to accommodate 
within our HQSS scheme. 
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